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Abstract. The ionization fraction is an important factor in determining the chemical and 
physical evolution of star forming regions. In the dense, dark starless cores of such objects, 
the ionization rate is dominated by cosmic rays; it is therefore possible to use simple analytic 
estimators, based on the relative abundances of different molecular tracers, to determine the 
cosmic ray ionization rate. This paper uses a simple model to investigate the accuracy of two 
well-known estimators in dynamically evolving molecular clouds. It is found that, although 
the analytical formulae based on the abundances of Hj ,H2,CO,0,H20 and HCO^) give a 
reasonably accurate measure of the cosmic ray ionization rate in static, quiescent clouds, 
significant discrepancies occur in rapidly evolving (collapsing) clouds. As recent evidence 
suggests that molecular clouds may consist of complex, dynamically evolving sub-structure, 
we conclude that simple abundance ratios do not provide reliable estimates of the cosmic 
ray ionization rate in dynamically active regions. 
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1. INTRODUCTION 

In this paper we examine the accuracy of the methods used to constrain the rate of cosmic -ray 
induced ionization in molecular clouds. The level of ionization within a dark cloud may deter- 
mine whether or not the cloud is stable against gravitational collapse leading to the formation 
of proto-stellar objects. This is particularly true if magnetic fields permeate the cloud - in which 
case the process is controlled by the coupling between the ions (tied to the magnetic field lines) 
and the neutral species. This applies both to quasi-static contraction (via ambipolar diffusion) 
and turbulent dissipation. Ruffle et al. (I1998> . investigated these processes and showed that the 
stability of low extinction (A^ < 3) clumps in the Rosette Molecular Cloud could be explained 
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in terms of photochemically driven high ionization levels, which inhibit the damping of (MHD) 
turbulent support. 

Cosmic rays also play a very important role in the physics of the interstellar medium - for 
example, they are an important source of heating in dark molecular clouds (e.g. Goldsmith & 
Langer [19781 1. These dark clouds are highly opaque to optical and ultraviolet (ionizing) inter- 
stellar radiation, yet they are known to contain a wide variety of atomic and molecular ions. 
Although dust grain surface chemistry is extremely important in dark and diffuse clouds (most 
significantly in the synthesis of H2) the chemistry within dark clouds is driven by gas-phase re- 
actions and the most significant of these (particularly for oxygen and nitrogen-bearing species) 
are ion-neutral reactions . The source of this ionization is believed to be cosmic rays. In the case 
of diffuse media the initiating reaction is 

H + crp — > + e 

(where crp stands for 'cosmic ray particle') 

whilst in dense molecular clouds approximately 97% of cosmic ray-H2 impacts lead to the 
formation of (Hartquist & Williams 1995) via the reaction: 



Subsequently, the then drives ion-molecule chemistries (see below). The rate of ionization and 
electron abundance within these clouds therefore has enormous implications for the efficiency of 
the molecular chemistry and the overall physical stability of many dark clouds. 

As a result of these reactions, nearly every cosmic -ray ionization creates an ion and this 
reacts with an abundant metal, such as oxygen or nitrogen (see below). The rate coefficient for 
reaction Q, ^, is therefore a key parameter in models of dark cloud chemistry; the relevant 
timescale for a cosmic-ray ionization induced chemistry to reach equilibrium is of the order 
of the fractional abundance of the metals (approximately 10"'*) divided by A value of f = 
1.3 X 10^'^s"' is often quoted in the literature, although in most cases only an order of magnitude 
accuracy is appropriate. In some sources, much higher values of ( have been inferred (e.g. at least 
^ = 8 X 10"'s ' in o Per, van Dishoeck & Black[l986). Generally, studies typically indicate that 
^ - 10"'^ - 10^'^s"' per hydrogen atom. 

However, if we assume a 'standard' value of ^ ~ 10^'^s"', the chemical timescale is ap- 
proximately 10^ years. Usually, this assumption of equiUbrium is adopted when using molecular 
abundances to infer the value of but here we note that the free-fall collapse timescale of a gas 
cloud is 3.4 X lO^n"'^ years (Spitzer [l978> so, for typical dark cloud densities of 10^ - lO'^cm"^, 
a dynamical timescale of 3.4 x 10^-10^ years is inferred. ( is usually estimated from observa- 
tions of the abundances of molecular species interpreted through simple chemical models. Such 



H2 + crp 




(1) 



which is followed by the rapid reaction: 



H+ -H H2 ^ -I- H 



(2) 
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methods depend on the assumption that chemical equilibrium has already been reached. It is 
not obvious, therefore, that estimates of the cosmic ray ionization rate in regions which may be 
rapidly evolving are correct. 

The rates of ionization have been clearly established for diffuse clouds, largely through ob- 
servations of OH and HD (see below). In dark clouds the situation is not so clear The standard 
assumption is that the same value of ^ is appUcable in dark clouds as in diffuse clouds. However, 
even in diffuse clouds, there is now believed to be very considerable variation in f from source 
to source. Thus, in a careful study based on H^ observations in ^ Persei, McCall et al. (I2003> . 
deduce that ^ is as high as 1.2 x 10^'^s"'. However, Le Petit et al. (2004) find that their model 
reproduces observed abundances in the same source using a value of ^ = 2.5 x 10"'^. In addition, 
as dark clouds are very much denser and also are permeated by ordered magnetic fields there is 
no a priori reason to assume that f is the same in diffuse and dark clouds. Moreover, as indi- 
cated above, attempts to estimate ( in dark clouds are usually based on estimators that make an 
assumption of chemical equilibrium. 

Most matter in a molecular cloud is essentially sterile, in that stars only form from a small 
fraction of the mass of a cloud (Leisawitz et al. I1999> . Most of this mass exists in regions of 
low extinction and is effectively in a photon-dominated region (Evans n"999> . Thus, most diffuse 
clouds are either supported against gravitational collapse, or else are transient objects that are 
dispersed before collapse and star-formation can occur This implies that dark, dense, clouds 
form from diffuse clouds and may retain some degree of the chemical identity of their precursors. 

In this paper we assess the implications of the assumption of chemical equilibrium on the 
inferred cosmic ray ionization rate in molecular clouds which are evolving on timescales that are 
short compared to the time required for chemical equilibrium to be established. 

2. ESTIMATES OF THE COSMIC RAY IONIZATION RATE, ^ 

The usual method of estimating the cosmic ray ionization rate is to employ a semi-analytical, 
semi-empirical technique. Thus what is, in reality, a complex multi-channel chemistry is sim- 
plified to to a small number of primary loss and formation reactions. The absolute abundances, 
or abundance ratios, of key tracers are then interpreted through simple analytical expressions. In 
practise, whilst many of the chemical abundances in these expressions are directly measureable 
(such as those for CO, HCO^ and Hp, others are unobservable (e.g. H2O) or else observational 
data does not exist. Consequently, the values used as input to the formulae are a combination 
of observed quantities and predictions from (often complex) chemical models. However, both 
the derivation of the analytical expressions and the modelled abundances used as input almost 
invariably depend on the assumption of chemical equilibrium. Therefore, for the reasons stated 
above, it is possible that they may not yield accurate values for ^. 

Lepp, Dalgarno & Sternberg ! 1 987l used the relationship between the abundances of (observ- 
able) OH and H3 and the equilibrium relationship between n(H3 ) and ^ (see below) to deduce 
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constraints on ^ in dense, dark interstellar clouds. This analysis was based on the assumption that 
the photodissociation of H2O by cosmic -ray induced ultraviolet photons is a major source of OH 
in dark clouds. 

In contrast, the abundances of carbon-bearing species in diffuse clouds are not sensitive to ^, 
and the the best estimators of ^ come from observations and models of the OH and HD abun- 
dances (e.g. van Dishoeck & Black 119861 Federman et al. I1996t . The entry into the oxygen 
chemistry is somewhat different than in dense molecular environments. The source of O^ is cos- 
mic ray ionization of hydrogen atoms and resonant charge exchange with oxygen atoms; 

H+ -H O ^ 0+ -i- H 

This is followed by reactions with H2 and dissociative recombination 
0+ -H H2, e" . . . ^ OH, H2O 
whilst the HD formation route is 
H+ -H D ^ H + D+ 
D+ H2 ^ HD H- H"^ 

Thus the abundances of both OH and HD are directly proportional to ^. Studies using this re- 
action network yield a value of ^ ~ 10"'^ - 10^'^s"^ for the rate of cosmic ray ionization of H 
atoms in diffuse clouds (Federman et al. ll996> - the spread of values being partly attributable to 
different assumptions being made and partly due to genuine variations, which may derive from 
uncertainties in the temperatures and the assumed radiation fields. 

If localised heating (e.g. through the action of shocks) is present, then the reaction 

O H- H2 ^ OH -I- H 

can occur, in which case these methods may overestimate the value of ^. 

An alternative estimator is based on the H^ abundances. In diffuse clouds, H3 is formed by 
reactions^andlJland lost by dissociative recombination: 

HJ -H e ^ H2 H- H 

although this approach can be undermined by uncertainties in the value of the dissociative re- 
combination rate co-efficient and the electron density. 

In diffuse clouds the main loss channel of simple molecular species is photodissociation 
and the implied timescales are of the order of 300 years, significantly less than the dynamical 
timescales of diffuse clouds. In addition the chemistry and rate-coefficients are are well-defined. 
For these reasons, the inferred ionization rates are believed to be reasonably robust (Hartquist 
and Williams II995I 1. However, McCall et al. 2003 present a direct detection of H:| in diffuse 
clouds towards ^ Per, and infer a correspondingly high cosmic ray ionization rate, 40 times 
larger than the 'standard' value. However, Le Petit et al. 2004, noting that such a high ionization 
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rate may be inconsistent with other molecular and atomic observations (especially those of OH 
and HS), explain the observed abundance of via a two-phase model. This example illustrates 
the difficulty of seperating the ionization rate from effects of structure and, as we discuss below, 
dynamical evolution). 

In dark clouds, the reactive species provides important entries into the ion-neutral chem- 
istry. H3 is usually more easily detected in its deuterated form, HiD^, to which it is linked via; 

-I- HD <^ H2D+ + H2 

This reaction also leads to the deuteration of molecular ions, such as DCO^, which can be used 
to constrain the fractional ionization (Caselli et al. 1998). 
The oxygen chemistry follows a simple route: 

-H O ^ OH+ + H2 (3) 
OH+ -H H2 ^ H2O+ + H (4) 
HzO^ +U2^ H3O+ + H (5) 
H3O+ -H e ^ OH -H 2H, H2O -I- H, OH -H H2 (6) 

Other major loss routes for H^ are reaction with water; 
Hj + H2O H3O+ -I- H2 (7) 
followed by dissociative recombination and reaction with CO, leading to HCO* formation; 
HJ -H CO ^ HCO+ (8) 

The HCO^ is lost through dissociative recombination; 
HCO+ + e ^ CO -H H (9) 
reaction with water; 

HCO+ -H H2O ^ H3O+ -t- CO (10) 
and charge transfer with low ionization potential elements; 

HCO+ +Na,Mg — ^ Na+,Mg+ + H + CO (11) 

This last reaction is particularly important in determining the ionization level of dark clouds. 
In metal poor dark clouds, HCO^ can be the dominant ion. However, the radiative recombination 
rates for elemental ions, such as Na^ and Mg^, are very slow so in undepeleted regions, reaction 
dl 11 1 can effectively result in a higher, stable, ionization level. 

WilHams et al. J 19981 1 used observations of C'^O, H'^CO+ and DCO+ towards 23 low-mass 
cores and a steady-state chemical model to deduce a mean value of ^ = 5 x 10 '^s '. The 
[DCO^]/[HCO^] ratio, which is relatively insensitive to ^ was used to deduce the ionization 
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fraction. In low mass cores, the ionization can be enhanced as a result of the presence of newly- 
formed low mass stars (which are X-ray sources) and the penetration of external UV, if the 
clouds are inhomogeneous. Bergin et al. ( I1999I I also used obervations of [DCO^]/[HCO^] to 
infer the fractional ionization within Orion massive dark cores, where cosmic rays are believed 
to be the dominant source of ionization. Both studies initially considered variations of ^ within 
a limited range (1 - 15 x 10 '^s '). This was constrained by the condition of thermal balance 
at the temperatures and densities of the sources being considered, recognising that cosmic rays 
are the dominant source of heating in starless cores. They subsequently adopted a value of ^ = 
5 X 10 '^s ' for all their modelled sources. 

However, these techniques are complicated by the fact that the fractional ionization and ^ 
are only partly independent of eachother, depending also on other free parameters such as the 
oxygen and carbon depletion, the abundance of the low ionization potential metals and the den- 
sity. Caselli et al. il998> performed a detailed study of the ionization fraction, X(e ), and ^ in 
dense cores. They compared simple analytical expressions for these quantities (dependent on 
the two abundance ratios; [DCO+]/[HCO+] and [HCO+]/[C'**0] and the free parameters listed 
above) with the predictions of a more elaborate chemical model, evolved to steady-state, and 
found reasonably good agreement, particularly for ^ < 10^'^s"'. In addition, they also noted 
that the abundance ratios do show considerable time-dependence on timescales < lO^years. As 
the sources that they have modelled are all quiescent cores (from Butner et al. I1995> they con- 
clude that they are older than lO^years and that chemical equilibrium can be assumed. However, 
the chemical model that they used assumed that the physical conditions are fixed. In a dynami- 
cally evolving system, the chemical timescales may be longer and the initial conditions less well 
constrained. 

van der Tak and van Dishoeck J2000> also modelled the HCO^ and H^ abundances (in steady 
state) on the assumption that the H2O and O2 abundances are low (due to freeze out) unless 
T>100K, when mantle sublimation occurs. They found that the value of ^ inferred from the H^ 
is systematically higher than that inferred from the HCO^ abundances. Noting that HCO-H traces 
dense gas, they speculate that the main cause of this discrepancy is the presence of intervening 
clouds along the line of sight, but again we postulate that a significant cause may be the non- 
equilibrium character of the chemistry. 

As described above, in contrast to diffuse clouds, the relevant timescale in dark clouds is 
approximately 10^ years, which may be significantly longer than the dynamical timescale in 
some scenarios. Moreover, there is increasing evidence to suggest that many of the structures 
in the interstellar medium may be more dynamically active than had been previously believed. 
Thus, for example, Ebnegreen (I1999> and Hartman et al. (1200 U have suggested that many of the 
translucent clumps detected in CO are transient. This has prompted plausible models of (MHD) 
transient structure formation (e.g. Falle & Hartauist l2002> . These models explain the presence 
of small, non-gravitationally bound clumps, which are found to have lifetimes of ~ 10^ years. 
This timescale is significant as it is less than the time taken for the complete conversion of H 
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to H2 and for chemical equilibrium to be established. The chemical implications of an initially 
H-rich chemistry were investigated by Rawlings et al. (I2002> who identified the presence of 
strong chemical signatures in such circumstances. In a cloud of density ~ 10"*cm"^ the dynami- 
cal (free-fall), chemical and freeze-out timescales are all comparable (~fewxlO^ years), so that 
the temporal and spatial variations of molecular abundances are expected to be sensitive to the 
boundary conditions (Rawlings l2003> . 

We suggest, therefore, that in dynamically evolving regions, the origin of the ionization can 
be ambigious. Specifically, a significant fraction of the ionization that is inferred for dark clouds 
may be a 'fossilised' remnant of diffuse cloud conditions. In other words the dark cloud car- 
ries chemical signatures of its more diffuse nascent state and has not yet established dark cloud 
chemical equilibrium. If this is true, then the inferred ionization rates for these objects may be 
incorrect. 



3. CHEMICAL INDICATORS OF ^ IN DARK CLOUDS 

The chemical networks within dark molecular clouds are well-known and, in their early stages, 
follow a reasonably linear chain of reactions. Without recourse to complex chemical models, 
it is possible to identify the key formation and destruction channels for those species which 
are sensitive to the the value of ^. This simplification allows simple analytical expressions which 
allow the determination of ^ from a knowledge of the abundances of the relevant chemical species 
to be defined. We investigate the validity of two well-known analytical expressions for f ; 

1. ^1, based on the abundance of H3 (see e.g. Lepp, Dalgarno and Sternberg ll987> . H3 is pri- 
marily formed by cosmic ray ionization Q and rapid reaction with H2 0. In dense environ- 
ments it is mainly lost through reactions with O, H2O and CO ((jSjl, 0, and (jSjl). Thus, in 
equilibrium: 

^1 = + 1^(0) + I^{H20)] (12) 

2. (2, based on the HCO^ abundance (see e.g. van der Tak & van Dishoeck l2000> . HCO^ is 
formed by reaction and primarily lost through dissociative recombination and reaction 
with water J10> . Thus, after substituting for ^(Hj ) from equation il2i we obtain; 

_ n(HCO^) [ifjf(CO) + 1^(0) + l^(H20)] 
l^(C0)X(H2) 

x[l^(e-) + l^{H20)] (13) 

In these expressions, the subscripts of the rate coefficients (fe,) refer to the reaction numbers given 
in the previous section, and the fractional abundance of a species, /, is defined to be the ratio of 
its number density to the total hydrogen nucleon density; 
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Species 



Abundance/X {H) 



He 



0.1 



c 



1.87 X lO-t 



N 



1.15 X lO-t 



O 



6.74 X 10-' 



S 



8.00 X 10 



Na 



2.00 X 10-" 



Table 1. Elemental abundances by number, relative to hydrogen 
4. THE MODEL 

The accuracy of both estimators for ^ was investigated using a one-point model of the time- 
dependent chemistry. The model treats protostellar clouds as being of uniform density and tem- 
perature, both of which change with time in three distinct phases: 

1. Phase I: Diffuse cloud conditions (no = lOOcm-^, T=100 K, A,,o=0.5). The system is allowed 
to evolve until no abundance changes by more than 0. 1 % in an individual time step. At this 
point we deem equilibrium to have been reached. 

2. Phase II: Spherically symmetric, homogeneous, free-fall collapse to a pre-specified terminal 
density (n„Mv - lO^cm"^). Initially, the collapse is assumed to occur isobarically, so that the 
temperature is inversely proportional to the number density, subject to the constraint that it 
does not fall below some minimum value (Tdark - lOK). As the collapse is assumed to be 
both spherically symmetric and homogeneous, the extinction A,, is given by 



3. Phase III: Relaxation to chemical equilibrium in dark cloud conditions (constant density and 
temperature, Tdark-^^ K) 

The chemistry includes 87 gas-phase and 36 solid-state chemical species, composed of the 
elements H, He, C, N, O, S and Na and linked through 1265 chemical reactions. The elemental 
abundances are given in Tabled These values are broadly representative of cosmic abundances, 
with a depletion factor of 0.5x included for carbon. (Recent observations (e.g. Meyer, Cardelli & 
Sofia 1997) have suggested that nitrogen and oxygen are also depleted by 50% in diffuse clouds. 
However, we find that including this extra depletion results in only extremely minor changes to 
our results.) 

We have investigated two extremes of metallicity; high (as in Tabled and low (as in Table[l] 
but with the abundances of S and Na (the representative low ionization potential metal) reduced 
by a factor of O.Olx. These values are consistent with those used by Caselli et al. J 19981 1 and 
others. 

The rate coefficents, kt, for key reactions /=2,3,7,8,9 and 10 are taken from UMIST ratefiles 
(Millar et al. 119911 Il997> and are given in Table |2 It should be noted, however, that since our 
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2.08 X lO^' 


h 


8.0 X 10-'" 




5.9 X 10-' 


ks 


1.7 X 10-' 


h 


1.1 X 10-' X (300/r) 


kio 


2.5 X 10-' 



Table 2. Rate coefficients for key reactions (cm^s ') 

study is concerned with the comparison of analytical indicators with modelled abundances, both 
of which depend on the same rate coefficients, our results are not sensitive to the adopted values 

of ki. 

The passage ffirough regions of intermediate extinction is assumed to be rapid and so no 
attempt has been made to model the details of ffie photochemistry, other than ffirough the normal 
dependence on A,,. We do, however, assume that H2 is self-shielding against photodissociation 
at all times, whilst photoionization of C and photodissociation of CO can only occur in diffuse 
cloud (phase I) conditions. 

Other than H2 formation (for which we assume k^u^, = 1.6 x IQ-'^cm^s-', Rawlings et al 
I1992> . in ffie standard models we have only considered gas-phase chemistry. However, in some 
models we have also allowed gas-phase species to freeze out onto the surface of the grains once 
ffie extinction rises above some critical threshold value (A,,c,;,. ~ 3). As the sources are dense and 
cold we do no need to consider the effects of surface chemistry as the return of species contained 
in ice mantles to ffie gas-phase is unlikely to be efficient. 

5. Results 

5.1. Evolution with time 

The evolution of selected species for a typical run is shown in Figure [0 This figure clearly 
shows how the the chemistry starts from primarily atomic values appropriate for diffuse cloud 
conditions. As the cloud collapses and becomes denser and more opaque to photodissociating 
radiation, the chemistry adjusts to a more molecule-rich, lower ionization state characteristic of 
dark clouds. 

In order to assess the validity of the analytical estimators of ^, we calculate the logarithmic 
ratio of the rates inferred from the two analytical expressions (^1, ^2) to the 'true' value used in 
the chemical models i^rme)' 

Ri = logio(-^) (16) 

Vfctrue / 

and we plot this value as a function of time. Thus, if ^, were to give a perfectly accurate value of 
^ ffien Ri would be zero. 

Results are given in Figure 13 for several values of ^,rue (10-'^ IQ-'^ 10-'^ IQ-'^ and 
lO-'V). 
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The first thing to note about these results is that, at late times (in phase 111, once dynamical 
activity has ceased) the curves all converge to values that are very close to 0. This is rather reas- 
suring in the sense that (although the estimators are based on huge simplifications of the chemical 
networks) the expressions are remarkably accurate indicators of the cosmic ray ionization rate 
for a wide range of ( when chemical equilibrium and dynamical steady state can be assumed. 

The second point to note is that whilst there are clearly very large discrepancies during the 
earUer (collapse) phase, the expressions give accurate results almost as soon as the collapse phase 
is halted, with no notable relaxation period. 

However, during the collapse (phase 11) we can see that significant discrepancies are present; 
(\ consistently underestimates (,rue, but for the lower ionization rates, this is only by a factor 
of ~ 2 - 3. However, if f„„(. were as high as 10"'^ - 10"'''s"', then the factor is much larger 
~ 10 - 10^. ^2, which is derived incorporating additional reactions in the same network used to 
derive ^i, is extremely sensitive to the abundances of species such as H2O, HCO^, and especially 
CO (See equation [O}. These abundances change dramatically during the collapse phase. The 
complex structure visible in Figure|3b) is partly attributable to this chemical complexity and is 
partly a result of inaccuracies in the numerical integration. However, as with (,\, the magnitudes 
of the discrepancies are largest in the ealiest stages of collapse. But, in contrast to ^1, R2 spans 
a smaller range (~ -1 to 1) for the values of ^tme that were investigated. Also we can see that 
the smaller values of ^,n,f (< 10^'^s"') are overestimated hy (2, whilst the larger values of ^„.„e 
are underestimated. From these results it would seem that (2 would always be > 10^'^s"', even 
though ^rrue may be somewhat smaller. 

5.2. Sensitivity to the free parameters 

We have investigated the sensitivity of our results to the various free parameters in our model. 
These include the terminal density (n,nax), the terminal temperature (Tdark), the choice of initial 
conditions corresponding to diffuse cloud equilibrium or atomic abundances, whether or not 
freeze-out processes are included, the metallicity, and the rate of collapse in phase 11. 

In models where the final density (n„,„x) was varied by up to two orders of magnitude, ^1 was 
slightly closer to zero for the smaller values of (,rue, whilst R2 was found to be quite insensitive 
to nf}^Qx. 

A range of post-collapse temperatures {Jdark) between 5 and 20 K were also investigated for 
a test case with fmif - 10^'^s"'. The effects on R\ were small (~ 20%) and are restricted to the 
early stages of collapse. In these early stages R\ is larger (closer to zero) whereas the opposite 
effect is noted for R2. The dependence on Tdark is even more marginal and is also restricted to 
early times, but here it is the lower temperatures that result in R2 being closer to zero. 

It has been shown (Rawlings et al. I2002> that the time-dependence of the chemistry may 
be sensitive to the assumed initial conditions - and in particular the initial H:H2 ratio. In our 
standard model we assume that the gas is initially characterized by diffuse cloud conditions, with 
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all hydrogen in molecular form. To investigate the possible effects of differing initial conditions 
on the validity of the ^ estimators we have also performed calculations where the gas is assumed 
to be intially atomic (H, He, C+, N, O, S+ and Na+) with no H2. 

We find that R2 is hardly affected by the choice of initial conditions, but some differences are 
noted for R\. Results are shown in Figure|3] For the larger values of (,rue, (i is larger by a factor 
of ~ lOx - bringing it closer to zero, in the early stages of phase II. For smaller values of (true 
(< 10 '^s^') the differences are rather marginal. 

Estimates of f and the fractional ionization are closely linked, and the latter will be coupled 
to the assumption that are made concerning the metallicity of the gas. Therefore, in addition to 
our standard abundances we have considered a "low metallicity" model, in which the abundances 
of the two low ionization potential metals (S and Na) are reduced by a factor of 100 to 8 x 10"^ 
and 2 x 10"*^ respectively. 

The results show a remarkable insensitivity to the metallicity and, as noted above, only a 
relatively small dependence on the initial conditions. This helps confirm the hypothesis investi- 
gated by this study; that the f estimators used are more sensitive to the dynamical status of the 
observed sources than any other individual chemical or physical parameter 

To simulate accelerated collapse (Lintott et al. I2005> we artificially increased the rate of 
change of the density by a factor of two. This also had very minimal affects which further sup- 
ports the idea that the observed discrepancies originate from near-instantaneous dynamical ef- 
fects rather than the existence of a longer-term chemical hysterisis. Other factors investigated 
were the extinction, and the temperature of the cloud prior to collapse. Neither had a significant 
effect. 

The result obtained from these parameter variations is therefore that R\ and R2 are relatively 
insensitive to the various free parameters and are primarily determined by the dynamical activity 
within the cloud. 

In our standard model we did not include any form of gas-grain interaction, but we have 
also investigated the effects of freeze onto the surface of the dust grains by allowing gas-phase 
species to freeze-out once A,, >3.0 mag. - similar to the critical value of the extinction required 
for the presence of ice water bands in Taurus (Whittet et al. 1200 U . This means that freeze- 
out commences during the collapse phase (II). No desorption mechanisms have been included. 
Results are shown for both Ri and R2 in Figure|3 

Two clear conclusions can be drawn from this figure: Firstly, if gas-grain processes are oper- 
ating - as is likely in the denser regions of the interstellar medium - then the analytical estimators 
always yield inaccurate values for ^, and secondly, the origins of the discrepancies (deviations 
of Ri from 0) are separable. In the collapse phase the differences have a dynamical origin, as 
in the gas-phase models, whilst in the static, post-collapse phase the differences are entirely a 
consequence of the freezing out of molecules onto dust grain surfaces. 
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6. Discussion and Conclusions 

Two related methods of estimating the true ionization rate, (,r„e, have been tested for a wide range 
of parameters. These were based on chemical networks centered on reactions of H3 ((2 using an 
extended version of the network used by ^1). In those cases where gas-grain interactions were 
ignored, we found that chemical equilibrium is rapidly re-established after the end of collapse, 
and both analytical expressions give accurate measures of the true cosmic ray ionization rate 
in the absence of dynamical activity. If freeze-out processes are important, then the analytical 
expressions are inaccurate at all times. 

There are two components to the modelled discrepancy between the analytical expressions for 
^1, ^2 and the true cosmic ray ionization rate ^tme- Firstly, the fact that the analytical expressions 
(^1 and ^2) are simplifications of a rather complex chemistry and only make allowance for the 
dominant chemical pathways, and secondly deviations from the assumption of equilibrium. Of 
course it is the latter that is the subject of this paper but, as noted in the previous section, it is 
reassuring to see that - at late times - the simple analytical expressions are reasonably accurate 
and the contribution from the former is relatively small. However, this is not always the case and 
has been investigated in some detail by Caselli et al. J1998t (see their Figure 9) who concluded 
that, even in an equilibrium approximation, the analytical expressions may overestimate ^ by a 
factor of ~ 2 - 3 for the cases where ( > 10^'^s"' . 

In our models we (obviously) use the same set of rate co-efficients in the analytical expres- 
sions and the full time-dependent models. Therefore, as the results are presented as the ratios 
of these quantities, they are reasonably robust and insensitive to the values of the adopted rate 
co-efficients. 

In all of the cases that we investigated, however, we found that the analytical expressions are 
very inaccurate in rapidly evolving, dynamically active regions. Discrepancies of several orders 
of magnitude are present, especially for the largest values of (true- Such conditions may, indeed, 
be appropriate in many dark clouds that have previously been thought of as quiescent entities. In 
support of this hypothesis, it should be noted that there is an increasing body of evidence that 
suggest that dark molecular clouds are not monolithic entities, but consist of ensembles of clumps 
and sub-structure that are not resolved with single-dish molecular Une observations. For example, 
populations of distinct sub-cores have been seen in TMC-1 Core D (Peng et al. ll998> and L673 
(Morata et al. I2003t . Moreover, most of these cores are not self-gravitating, but are probably 
transient entities or varying size/mass and dynamical status. This picture has been supported by 
recent chemical models (Garrod, Williams & Rawlings I2005t which identify the sub-cores as 
objects that grow and decay on timescales of ~ lO*" years. 

Some sensitivity to physical parameters was found, in particular to the density reached after 
collapse. However, these effects were much smaller that the differences observed in the same 
core at different stages of collapse. It can therefore be seen that these methods cannot accurately 
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determine the cosmic ray ionization rate unless collapse can be definitely shown to have been 
completed. This forms a challenging observational requirement. 

Geballe et al. 2003 use the observed abundance in several dark clouds to constrain the 
cosmic ray ionization rate. Oka et al. 2005 compare these measurements of cosmic ray ionization 
to their measurements of ^ in diffuse clouds. However, the results presented in this paper suggest 
that such a comparision is invalid unless the dark clouds being studied were sufficiently old for 
the dynamical effects revealed by our models to be unimportant. 

In their study of massive star-forming regions, van der Tak and van Dishoeck (I2000t found 
that the cosmic ray ionization rate as deduced from H3 observations was greater than that as 
deduced from HCO^ observations. That is to say, models using ^ derived from H3 observations 
resulted in a significant over-production of HCO^. They speculated that, whilst the HCO^ mea- 
sure probes only the dense molecular gas, contributions from foreground layers (e.g. in photon 
dominated regions or translucent clouds) may be responsible for the discrepancies. This is a 
valid explanation, but we note that dynamical activity may help explain the discrepancies; the 
differences between the values of ^ inferred from the H3 and the HCO^ observations (effectively 
quantified by Ri - R2) depends on the 'true' value of ^ and the dynamical status of the observed 
source. The more dynamically active the source, the greater the discrepancy. 

In conclusion, it is possible that, if a reasonably full dataset of observations of key molecular 
tracers is available for any given source, the absolute values of ^1 and (2 and the difference be- 
tween .Ri and R2 (and perhaps other analytical expressions for f ) could be used to simultaneously 
constrain the true ionization rate and the dynamical status of the source. If it turns out that very 
similar values are obtained for ^1 and ^2 then that is an indication that the observed cores are 
chemically 'old' and dynamically quiescent. 
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